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PREPACE 


Part  of  the  Project  RAND  research  program  consists 
of  basic  supporting  studies.  The  mathematical  reseaioh 
presented  here  Is  concerned  with  an  analytic  technique 
for  approximating  a  given  function  by  a  svun  of 
exponentials.  The  results  that  are  obtained  have 
applications  In  control  theory,  clj?cult  synthesis,  and 
numerical  analysis. 
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SUMMAKY 


;Ve  consider  a  simple  analytic  technique  for  flndinc 
linear  differential  equations  which  are  approximately 
satisfied  by  a  clven  function,  and  thus  obtain  approxi¬ 
mations  to  a  civen  function  by  a  sum  of  exponentials. 
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A  NOTE  ON  DIFFERENTIAL  APPROXIMATION  AND 
ORTHOGONAL  POLYNOMIALS 


1 .  INTRODUCTION 

The  problem  of  obtaining  an  exponential  polynomial 
N  X.t 

of  the  form  2  b.e  vrtilch  closely  appixjxiraates  a 
1=1  ^ 

given  function  f(t)  In  an  interval  a  ^  t  <  b  is  a 
problem  of  some  difficulty  if  we  allow  both  the 
coefficients  and  the  exponents  to  be  unknoms.  Either 
of  the  criteria  of  fit, 

N  X.t 

(1.1)  max  lf(t)  -  Z  be  ^  |, 

a^t^  1=1  ^ 


or 

fb  N  X.t  2 

(1.2)  /  |f(t)  -  2  b.e  ^  1  dt 

1=1  ^ 


lead  to  difficulties;  see  Lanczos  [  l]  . 

In  this  note,  we  ;d.3h  to  consider  a  different  way 

of  measuring  the  closeness  of  f(t)  to  a  sum  of 

exponentials.  If  r(t)  v:ere  a  function  of  the  fonri 
N  X  t 

2  b.e  ^  ,  it  would  satisfy  a  linear  differential 
i=l  ^ 

equation  of  the  form 


(1.3)  f^^^^  +  +  •••  +  Cj^  =  0. 


Hence,  let  us  attempt  to  determine  real  coefficients  c 
which  minimize  the  integral 
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(1.1)  /  °°  (f  <«>  +  *  .  .  .  +  o„f  )2dt. 

—00 

We  call  this  differential  approximation. 

Results  of  the  type  we  obtain  have  applications  In 
the  fields  of  control  theory,  circuit  synthesis,  and 
numerical  analysis. 


2.  FOURIER  TRANSFORMS 

It  Is  clear  that  we  must  Impose  some  conditions  on 
the  function  f(t)  In  order  to  pose  the  problem. 

Assume  then  that  e  L^(— oo,oo)  for  k  =  0,1,2,...,N, 

Then  If 


(2.1)  s(s)  =  r  f(t) 


f(t)  =  ^ 


— oo 


oo 


e  dt , 


c(s)e  ^®^ds. 


— oo 


we  have 


(2.2) 


.(k) 


(t)  =.  J  (-is)^G(s) 


-Ist , 
e  ds . 


— oo 


Using  the  Plancherel— Parseval  fonuula,  we  have 


(2.3)  f  (f^^^  + 


.(N-1) 


+  c 


dt 


— oo 
oo 


lc(s)  l^lcj^  -  Isc.^^  +  (-is) 


2  ,  |2  , 
Cjj_2  +  •  •  •  I  ds 


— oo 


/°°iG(s)|2[( 


'I'l-l  ''N-3 


— oo 
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3.  ORTHOGONAL  POLTOOMIALS  AI'TD  MINIMIZATION 

It  Is  clear  that  a  chance  of  variable  reduces  the 
problem  to  that  of  finding  the  orthogonal  polynomials 
associated  respectively  I'ri.th  the  weight  functions 
|c(s)|  and  s  lc{s)|  .  Since  these  orthogonal  poly¬ 
nomials  can  be  constructed  In  a  systematic  fashion,  we 
have  a  simple  v;ay  of  obtaining  the  coefficients  c^, 
and  further  measures  of  the  asymptotic  behavior  as 
N  -»  00  ;  see  Szego  1 2)  . 

A  direct  treatment  of  the  question  of  minimizing 
the  quadratic  form  of  (1.4)  would  not  be  simple 
computationally  for  large  N  (although  a  Schmidt 
orthogonallzation  could  be  used),  and  would  not  readily 
furnish  asymptotic  behavior.  The  technique  presented 
above  is  most  useful  In  connection  with  the  treatment 
of  vai’lous  classical  functions  where  s(s)  has  a  simple 
tmalytlc  form. 

Similar  results  can  be  obtained  for  the  case  where 
the  interval  is  finite  or  semi— Infinite,  but  not  of  the 
same  simplicity. 

4 .  STABILITY 

Tlie  study  of  the  precise  connection  between  the 
solution  of  the  linear  differential  equation 

(4.1)  u^^^  +  +  •••  +  CjjU  =  0, 


and  the  function  f(t)  leads  to  a  stability  question 
vdilch  we  shall  study  elsevrtiere. 
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